On spontaneous photon emission in collapse models 

Stephen L. Adle 

Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, USA 

Angelo BassJll and Sandro Donad 

Department of Physics, University of Trieste, 
Strada Costiera 11, 34151 Trieste, Italy and 
Istituto Nazionale di Fisica Nucleare, Trieste Section, Via Valerio 2, 34127 Trieste, Italy 

We reanalyze the problem of spontaneous photon emission in collapse models. We show 
that the extra term found by Bassi and Diirr is present for non-white (colored) noise, but its 
coefficient is proportional to the zero frequency Fourier component of the noise. This leads 
one to suspect that the extra term is an artifact. When the calculation is repeated with the 
final electron in a wave packet and with the noise confined to a bounded region, the extra 
term vanishes in the limit of continuum state normalization. The result obtained by Fu and 
by Adler and Ramazanoglu from application of the Golden Rule is then recovered. 



I. INTRODUCTION 



In a previous series of articles [lH3j , the problem of the spontaneous emission of radiation from 
charged particles, as predicted by collapse models, was analyzed in detail. The interest in this 
kind of problem arises from the fact that it currently sets the strongest upper bound on these 

-r fil — rr- M: — tt — — 

of [3] has been done within the QMUPL model [7|]. In spite of the fact that the two models should 
give the same predictions in an appropriate limit [3j|, the formula derived in 131 turns out to be 

n q 

twice bigger than that of pj, |2J], in the case of radiation emitted from a free particle. While this 
factor of 2 difference is unimportant in the free particle case and for a white noise, as we will show 
here it gives rise to awkward terms in the case of a colored noise, leading one to suspect that the 
extra term is not physical. 

We wish to come back on this issue, in order to clarify some mathematical details regarding the 
derivation of the radiation formula. 
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II. THE CSL MODEL FOR CHARGED PARTICLES 

In the CSL model, the standard Schrodinger equation is modified by adding terms which cause 
the collapse of the wave function: 



d\ik 



±Hdt + ^ [ dx [M(x) - (M(x)) t ]dWi(x) _ ^_ f dx [M(x) - (M{x)) t } 2 dt 
h m J 2mg J 



<?(x) = , J ■ e- x2 /24 (3) 



(1) 

here H is the standard quantum Hamiltonian of the system and the other two terms induce the 
collapse. The mass mo is a reference mass, which is taken equal to that of a nucleon. The parameter 
7 is a positive coupling constant which sets the strength of the collapse process, while M(x) is a 
smeared mass density operator: 

M (x) = J2 m J N J ( x ) . N i W = / 9 (x - y) ^ (y) if>j (y) d s y, (2) 

3 

V'j (y - ) s )) (y> s ) being, respectively, the creation and annihilation operators of a particle of type 
j at the space point y. The smearing function <?(x) is taken equal to 

1 

where rc is the second new phenomenological constant of the model. Wt (x) is an ensemble of 
independent Wiener processes, one for each point in space, which are responsible for the random 
character of the evolution; the quantum average (M(x))t = (ipt\M (x)\ipt) is responsible for its 
nonlinear character. 

As shown e.g. in Q], the averaged density matrix evolution associated to Eq. (pQ) can also 
be derived from a standard Schrodinger equation with a random Hamiltonian. Such an equation 
does not lead to the state vector reduction, because it is linear; nevertheless, since they both 
reproduce the same noise averaged density matrix evolution, and since physical quantities like the 
photon emission rate can be computed from the noise averaged density matrix, the non-collapsing 
equation can equally well be employed to compute such quantities. The advantage of this second 
approach is that, being based on a linear (stochastic) Schrodinger equation, it is much easier from 
the computational point of view. In our case, the stochastic Hamiltonian is: 

H TOT = H-h^J2 1 ^ [ iV(y,t)Vj(y)^(y) d 3 y (4) 

3 

where: 

N(y,t) = j g(y-xMx)d 3 x, (5) 
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and £t(x) = dWt(x)/dt is a white noise field, with correlation function E[£ t (x)£ s (y)] = 5(t — 
s)S(x. — y). As such, iV(x,i) is a Gaussian noises field, with zero mean and correlation function: 

E[N( X ,t)N(y,s)] = 5(t-s)F(x-y), F(x) = —±—^l^c. (6) 

The purpose of this article is to reconsider the analysis of the emission of radiation from a free 
charged particle, previously discussed in the literature. Accordingly, in the following we will be 
interested only in one type of particle, so from now on we will drop the sum over j. 

The Hamiltonian H TOT can be written in terms of an Hamiltonian density H TO t- For the 
systems we are interested in studying, we can identify three terms in M TOT : 

^tot = H p + H r + Hi N t- (7) 

H P contains all terms involving the matter field, namely its kinetic term, possibly the interaction 
with an external potential V, and the interaction with the collapsing-noise: 

vn 

Up = — VV>* • VV> + Vip*ip - hyFj — Nip*ip. (8) 
2m mo 

Hn contains the kinetic term for the electromagnetic field: 

w »=K eoEi+ £>' (9) 

where Ej_ is the transverse part of the electric component and B is the magnetic component. 
Finally H INT contains the standard interaction between the quantized electromagnetic field and 
the non-relativistic Schrodinger field: 

Hint = t — V>* A ■ Vt/> + — A 2 ^>- (10) 
m 2m 

The electromagnetic potential A(x, t) takes the form: 

A(x, t) = Y, « P [e P ,A a p e^-^ + e PjA a^e'^-^] , (11) 

p,A 

where a p = ^h/2eou> p L 3 and oj p = pc. We are quantizing fields in a cubical box of size L. 

To analyze the problem of the emission rate, we will use a perturbative approach. We identify 
the unperturbed Hamiltonian as that of the matter field (interaction with the noises excluded) plus 
the kinetic term of the electromagnetic field: 

Ho = ^Vf • VV + VW + H R , (12) 
2m 
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and we assume that its eigenstates and eigenvalues are known. In particular, we assume that the 
matter part %q is diagonalizable. The perturbed term then is: 

U x = i—ip*A-V?p + — AVV>-V7 — NiP*i>. (13) 
m 2m rriQ 

Such a division of M TOT in %q + Hi is justified by the fact that the effects of spontaneous col- 
lapses driven by the noise field are very small at microscopic scales. This is also true for the 
electromagnetic effects we are interested in computing. 

III. FEYNMAN RULES 

The Feynman diagrams for our model can be derived in a standard way, by means of the Dyson 
series and Wick theorem. They are: 

1. External lines (the symbol • denotes the generic space-time vertex (x, t)): 
k » . = u k (x)e-i Ekt • > k = u* k (x)eTi Ekt 

= JV(x,t) 

The functions u&(x) are the eigenstates of — ^V 2 + V, and E k is the associated eigenvalue: 



2m 

^ 2 + v 



itfc(x) = E k u k (pc). 



Since the noises field N is treated classically, there is no distinction between incoming and outgoing 
lines. 

2. Internal lines. The propagators for the matter field and for the photons, are: 

• > • = Fl2, • r \S\j* = P\2 ) 

12 12 

with 1 = (xi,ti), 2 = (x 2 ,i 2 ) and: 

F 12 = F(xi,ti;x 2 ,t 2 ) = e(t 2 -ti)^n fc (x 2 )<(x 1 )e-s^('2-*i) (14) 

pirn = p^ (xi)ti;X2)t2 ) = e(tl _ t2) ^ a 2^ Ae j.m e i[k.(x 1 -x a) -a, t ( tl - ta )] 

k,A 

+ « (t 2 - ii) ^«fc e k!Aek,Ae i[k ' {x2_Xl) ^ fc( * 2_<l)] - (15) 

k,A 

3. Vertices. There are three types of vertices, corresponding to the three terms in the interaction 
Hamiltonian T-L\\ 

— > ' > — m / \ m — * — ■— » — mo 
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In the first vertex, the derivative acts always on the incoming external line. In the second vertex, 
e 2 /m appears in place of e 2 /2m (as one would naively expect by inspecting at Eq. (I13p ) in order 
to take properly into account the multiplicity of the diagrams. The same rule applies also to the 
standard scalar QED (without the noise term). 
4. One has to integrate over space and time in all vertices 




Note that there is no factorial term 1/n! coming from the Dyson's series, because this is canceled 
by the multiplicity of the diagram 1 . Only diagrams containing double photon propagators, like: 




do not follow this rule. In such a case, one has to multiply by a factor 1/2 for each such loop. 



IV. PHOTON EMISSION PROBABILITY AT FIRST PERTURB ATI VE ORDER 

At first order in y^y and e, there are two different contributions to the process of photon emission, 
coming from the interaction of the free particle with the noise field: 



Y 2 1 V X 2 : \ 

According to the rules previously outlined, the contribution of the first diagram is: 

y 777, J y 777-0/ ^ Jt% J ti 

x[ dux f dx 2 n i (x 2 )e- ipXl e; A - [Vn fe (x 1 )]<(x 2 )n* / (x 1 )Ar(x 2 ,t 2 ), (16) 



while the contribution of the second diagram is: 



y 777 J y / r itj J ti 

x j rfxi jf 3 dx 2 n fe (x 1 )<(x 2 )e- i P X2 e; iA • [Vn,(x 2 )]n}(x 1 )iV(x 1 , (17) 



1 More precisely, a diagram containing n vertices has a factor 1/n! in front, coming from thew Dyson's expansion. 
However, there are n! such identical diagrams, differing only in the way the vertices are numbered. 
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Summing these two contributions, the transition amplitude Tfi becomes: 

( /| ^.'.e-'P'e;, . V |*) (fc| JV(x,t 2 ) |>) + (/| JV(*,tO |*) (*| e^e-**^ • V |i>] ,(18) 

where we have introduced the position operator x. It is convenient to rewrite the above expression 
in a more compact form. Since the correlation function © of the noise is a product of its time and 
space components, as far as the average values are concerned we can replace N(x,t) with ^ -/V(x), 
where & is a white noise in time, while N(x) is a Gaussian noise in space, with zero mean and 
correlator F(x. — y). We also introduce the following two operators: 

W = a p (i—) e-'P- x e p , A -V, Af = -h^j— iV(x). (19) 

\ m J mo 

The first operator refers to the radiative contribution (hence the symbol 1Z), the second one to 
the interaction with the noise (hence the symbol AT). Defining moreover the matrix elements 
lZ p ki = (k\1Z p \i) and Mki = (A;|jV|i), considering photons with linear polarization (e* A = e p ,\), we 
can write Eq. (|18p in the following way: 



e^^n p fk M hi + e^^MfkKl] . (20) 

This is the final expression of the first-order transition amplitude for a charged particle to emit a 
photon, as a consequence of the interaction with the noise field. The particle might be free, as we 
will consider in the next section, or interacting with an external potential. 

V. EMISSION RATE FOR A FREE PARTICLE 

In the case of a free charged particle, the initial and final states and the generic eigenstate of 
Hp are: 

Uj(x) = ; , Uf(x) = ; , Itfc(x) = . , (21) 



and we have chosen the reference frame where the particle is initially at rest. The corresponding 
eigenvalues are E\. = 7i 2 k 2 /2m, and similarly for Ei and Ef. The matrix elements for the radiative 
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part can now be easily computed: 



!< -fk 



(m p \k) 



L 3 



dx e 



-iqx 



L3 



a n I i — e 



-ipx 



e p,A 



„ikx 



Mi 



ki 



'Hi 



m L 3 7 



A/". 



(fc|7£ p |i> = 0, 

-h^/\— -L [ dx7V(x) 
m L s J 



-tk-x 



3 i(k-q)-x 



(22) 

(23) 
(24) 
(25) 



As we see, the contribution given by the second Feynman diagram is null. Therefore, in squaring 
Eq. ([20]) . taking the average with respect to the noise, we obtain the relatively simple expression: 



El 



T M 2 = ^EE^ nP fk m^ 



(26) 



k j 



x / dt\ I dt 



dt 3 

o 



t 3 



dt 4 e iatl e ibt2 e ict3 e idu 5(t 2 -t 4 ), 



io Jo 

where we have set U = and tt = t, and moreover we have defined the constants: 



a = - (E f + hw p - E k ) , b = -(Ek — Ei) , c . 



\{E f + hu p -Ej), d^ — iEj-Ei), 



n 



(27) 



We focus the attention on the temporal part. We have: 



o Jo 

t 



JO 



dt x I dt 2 I dt 3 I dt^ atl e ibt2 e %ct3 e idu 5{t 2 -U) 



[ dh [ dt 2 [ dt 3 [ dUe^e^e^e^Sih-t^eih-t^eits-U) 
Jo Jo Jo Jo 

dt 2 [ dh [ dt 3 e iatl e iib+d)t2 e ict3 

J to J to. 



1 

ca 



J(c+a)t 



1 - e i9t 



+ e 



i(g+c)t _ i p i{g+a)t _ i i _ i(g+c+a)t 

iat^ i ret c i 



i o + c 



i(g + a) i(g + c + a) 



(28) 



where we have defined: 



g = b + d = -(Ek-Ej) 



Because of the relation a + b + c + d = a + c + g = 0, Eq. ([28]) simplifies to: 



1 



ac 



e -igt _i e 
: + 



— + : t 



(29) 



(30) 



ig la ic 

We are now ready to replace the matrix elements 1Z P * an d 7£j fc with the explicit expressions ([22]) 
for the free particle. The indices k,j become vector indices k,j labeling the wave number, and the 
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constraints given by the deltas in the 1Z terms (see Eq. (|22p) suppress the two sums in Eq. ([26]) . 
This moreover implies: 

g = (31) 

a = " c = \ ( E f + ~ E i+p) = {p c -^~ ' (32) 
Then expression (|30j) for T simplifies to: 

2 

T = -3 [at - sin (at)] , (33) 
We now focus on the remaining part of Eq. (|26p : 

^EE ^ p fk n^ kl ] = ^1 (^) 2 ( £p ,a • q) 2 n^ +q) ^ Mi ], (34) 

where we have taken into account the constraints coming from the Kronecker delta in Eq. (I22p . 
The stochastic average gives: 

E[Ar* p+q)t Af (p+q)l ] = n 2 7 -1 ^ dxx ^ dx 2 e'CP+D-C^^ir ( Xl - x 2 ) . (35) 

We make the change of variable: x = xi — x 2 and y = xi + x 2 (the Jacobian is 1 /8) and we use 
the rule: 

r+f /■+§ 1 r+L f + (L-x) 



f+2 f+2 1 f+ L f + ^-x) 

/ dxi / dx 2 f(xi,x 2 ) = - / cte / dy [/ (x, y) + / (-z, y)] , (36) 

J-\ J-\ 2 JO J-(L-x) 

thus obtaining: 

= h 2 j(—) 2 17 / +L ^i l HL ~ Xi) dyi— r Lfei^i+e-i^ile-^o. (37) 
Vm y 8L b f = iJ J-(L- Xt ) V4Trr c l J 



The integral over y, gives: 



2 ^ J-(L-xi) L 
The second term vanishes in the large L limit, so we can ignore it. We are left with: 

m^^Mi] = ^(^) 2 p /_ + L L ^^)- x F(x). 



^ r+(L-Xi) 

d Vi = 1 - (38) 



(39) 



In the large L limit, the integral gives the Fourier transform of the correlation function F. Taking 
into account the form ([6]) of F, and collecting all pieces, we have: 



E\T f ^ = A(e p , A • q) V«^)H ^frO (40) 
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where: 

collects all constant terms. 

The emission rate T(p) can be computed from the transition probability, by differentiating over 
time, and by summing over the momentum q of the outgoing particle and the polarization A of 
the emitted photon, according to the formula: 

A 

Let us choose the axes so that p = (0, 0,p); in this way a, as given by Eq. (|32p . becomes a function 
only of p and q z . The sum over polarizations then gives ^a( 6 p,a " °l) 2 = Qx + Qy- All factors L 
cancel with each other, so we can take safely the limit L — > +oo. The sum over q then becomes a 
triple integral. The two integrals over q x and q y can be easily computed, being Gaussian, and we 
obtain: 



J^K£)(i)/^~ fe+p,vsi ^' (43) 

The above integral can be rewritten in the following way: 

dq z e-^^ 2 c l - C f at) =^J dze^ 1 = , (44) 

where we have defined the following new quantities: z = hp(q z +p)/m, D = pc + hp 2 /2m and 
j3 = mr c /hp. Since /3 ~ 1CP 13 s and D ~ pc ~ 10 19 s _1 for a non relativistic electron and 
for radiation in the KeV region, the Gaussian term in Eq. (|43p is vanishing small in the region 
where 1 — cos[(D — z)t]/(D — z) 2 is most appreciably different from zero. Around the origin, 
where the Gaussian is not negligible, the denominator varies slowly, and one can approximate 
1/(D — z) 2 ~ l/D 2 , and bring it out of the integral. What remains, apart the Gaussian term, is 
1 — cos[(D — z)t\. The second term oscillates very rapidly and gives a negligible contribution to 
the integral. Thus only the first term survives. When integrating moreover over all directions in 
which the photon can be emitted, the emission rate becomes: 

dT _ Xhe 2 

dp 2^ 2 £Q(?m\r 2 p'' 

r-i 

with A = 7/87r 3 / 2 7'e equal to the collapse rate first introduced in the GRW model [9J]. In the above 
expression, we have neglected the oscillating term, which averages to zero over typical experimental 
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times. The above result is expressed in SI units. The transformation to CGS units simply requires 
the replacement Eq — > 1/47T, in which case we obtain twice the results of 2j. 

The mathematical reason for such a difference lies in the type of approximations used to obtain 
the final formula. Going back to Eq. (|30p . in [2] the following approximation was made: 



1 

ac 



: + : + : t 



-• (46) 
ac 



ig la ic 

While this is legitimates in general, it gives problems in the free particle case. Here, as we have seen, 
g = 0, meaning that the oscillating term depending on g becomes linear in t. This contributions 
sums with the other linear term, giving the factor of 2 difference. We also note that the remaining 
two oscillating terms are mathematically important, though physically negligible. Since for a 
free particle a = — c, they reduce to a cosine, which makes sure that the integral in Eq. (|43p is 
convergent. Without it, the pole at the denominator would produce a divergence. 

VI. EMISSION RATE IN THE NON- WHITE NOISE CASE 

The factor of 2 difference, while not being significant in the white-noise case, has serious con- 
sequences in the colored- noise case. To see this, we now generalize Eq. (|45p to the case where the 
collapsing noise has a correlation function which is not white in time: 

E[N(x,t)N(y l8 )] = /(i-s)F(x-y). (47) 

We can start from Eq. (|26p for the average transition probability: 

E I T /*I 2 = i£E / dz n%M kl ^)n))Mp^) (48) 

x ( dti ( 1 dt 2 ( dt 3 ( 3 dt 4 e mtl e lbt2 e lct3 e idu f (t 2 - U) , 
Jo Jo Jo Jo 



where now / replaces the Dirac delta. The coefficients a,b,c and d are the same as in (127p . The 
only effect of the non-white noise is to modify the time dependent part of the transition probability, 
which we consider separately: 



T = [ dti j 1 dt 2 [ dt 3 [ ^ dt 4 e tatl e ibt2 e ict3 e idt4 f(t 2 -U). 
Jo Jo Jo Jo 



(49) 



This can be rewritten as follows: 

r-t ft ft ft 



T= f dt 2 f dh ! dt 4 f dt 3 e iatl e ibt2 e lct3 e idt4 f(t 2 -t 4 ) 
Jo Jt 2 Jo Jti 

= -—/"* dt 2 f dt A (e iat - e iat2 ) (e ict - e icU ) e ibt2 e idu f {t 2 - t A ) . (50) 
ac Jo Jo 
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t ft 

du I dve iau e iliv F{u-v) = I du I dv e* [(a+PX"*>)+(<*-PK«-v)] f ( u _ v ) . (51) 



There are four terms, which all have the following structure: 

du 

Jo Jo Jo Jo 

We perform the change of variable x = u — v and y = u + v. In these new variables, the integral 
changes as follows: 



t rt 

du j dvf(u,v) 
u Jo 



j /-t f2t—x 

2 J dx J dy if(~ x >y) + /O^y)] • 



(52) 



In our case, the integrating variables separate, and the integral over y can be easily performed, 
giving: 

*2t-x 



I 

J x 



dye 



e i("+/3)* 1 
4———sm-(a + P)(t-x). 
(a + p) 2 



(53) 



Then, taking into account that f(u — v) = f(x) = f(—x), the double integral / reduces to: 

e f(a+/3)i ft I 1 

J = 4- — — — - / dxf(x)sm-(a + P)(t-x)caB-(a-P)x. 
(a + p) Jo 2 2 

Going back to Eq. ([50]) . taking into account the relation a + b + c + d = 0, we can write: 



(54) 



4 I e~^ gt f l 
T = < / dx f(x)sm 



ac y g J 

& \at rt 



a Jo 

e \ct ft 



dx f(x) sin 
dx f(x) sin 



2^(* - x ) 



-ait — x) 
2 v ' 



-c(t — x) 



cos 



COS 



-(b-d)x 



1 

— ax 
2 





"1 


cos 


-ex 




2 



1 /■* 

H — / cix f(x)(t — x) cos [(a + 6)x] 

2 jo 



(55) 



Of course, in the white noise case f(x) = S(x), the above equation reduces to Eq. ([30]) . 2 In 
computing the matrix elements 7£?- and Ajj for the free particle, a further constraint comes from 
the Kronecker delta of Eq. (|22p . which implies a = —c and g = 0. Accordingly, the expression for 
T further simplifies to: 

2 ' 



, w fix f(x)(t — x) [cos(6x) + cos[(a + o)xl] 
a Uo 



4 fl \ I* 

— cos [-at] / dx/(x)sin 

a V 2 / Jo 



1 / 

-a(t-s) 



cos 1 -ax 



(56) 



Note that J Q da;5(a;)5(a;) = ^g(O), for a general function g(x), must be used in the reduction to the white noise 
case. 



12 



The next step, in computing the emission rate, is to compute the time derivative. Differentiating 
with respect to the upper limit of the integrals, produces terms proportional to f(t), which vanish 
in the large time limit, as we assume that the correlation function has a finite correlation time. 
The remaining terms coming from the second line produce oscillating terms, which average to zero. 
Thus, the only significant term, in the large time limit, is: 



dT 



1 



f(b) + f(a + b) 



dt t — oo a 2 L J 

where we have defined the Fourier transform of the correlation function: 



(57) 



7M = 2 



dt f(t) cos(wt) 



dtf(t)e 



iiut 



(58) 



Finally, in computing the integral over the final momentum of the particle (see Eq. (|43p ). we have 
have approximated the Gaussian term by a Dirac delta, meaning that we are imposing q ~ — p. 
This implies: 



pc 



hp 2 hp ■ q 

2m m 
/i(q + p) 2 



pc + 



2m 



pc 



2m 



0. 



(59) 
(60) 



Thus we have: 



dT 
dp 



NON-WHITE 



/(o) + f( P c) 



dT 

dp 



(61) 



The second term is the expected one: the probability of emitting a photon with momentum p is 
proportional to the weight of the Fourier component of the noise corresponding to the frequency 



pc. The first term instead is independent of the photon's momentum, and looks suspicious. 



Precisely this term, in the white-noise limit, is responsible for the factor of 2 difference, as one can 
easily check. In the remaining sections we analyze the origin of such an unexpected term. 



VII. COMPUTATION USING A GENERIC FINAL STATE FOR THE CHARGED 

PARTICLE 

At the end of Section [V] we have discussed that the factor of 2 difference arises because g as 
defined in Eq. ([29]) becomes 0, the reason being that the deltas in the 1Z terms (see Eq. ([22]) ) force 
Ek to be equal to Ej. One could then expect that by considering a generic final state for the 
outgoing particle — in place of the more artificial plane wave — such a constraint is removed. We 
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discuss such a possibility in this Section. So, instead of a final state for the particle with definite 
momentum, let us now take a normalized wave packet: 



1 e «(q+A)-x e ik-x 

n i (x) = - 7 =, u/(x) = ^ft(A) ^_ , u fc (x) = - 7 =, (62) 



A 

where (A) normalizes the wave function: 



The matrix elements (122I)-([25"D now become: 



(63) 



n% = (f \K p \ k) = £V (A) (/a |^| k) = a p (-^) J>* (A) [e p • (q + A)] £ k , q+A+p . (64) 

M ki = -hV\—^ [ (ixiV(x)e- ik - x (65) 
m L A J 

K p ki = (k\KP\i) = 0, (66) 

M fk = -hVX-^Yh(A) [ dxiV(xW( k - t i- A )- x (67) 

Since also in this case 7£?- = 0, the formula of E|Tjj| 2 is still given by Eq. (|26p . where the temporal 
part takes the same form as in (130 j) : 



e i t /*i 2 = ^EE n % n % E WM T > (es) 



with T given in Eq. (|30p . The two Kronecker deltas coming from Eq. (|64p set: k = q + A + p and 
j = q + A' + p. Accordingly, the coefficients a, c and g, defined in ([27]) and ([29]) . become: 

_ Cg/ + fr^p - E q+A+P ) _ (E f + frwp - Eq +A i+p) _ E q+A+P - -Eg+A'+p , , 

a ~ h > c - n ,9- n ■ (byj 

Moreover, we have: 

E[A^AT fci ] = n 2 7 f — V -L / rfx! / dx 2 e- l (P + ^ ( Xl - X2 )e- J ( A xi - A '- X2 )F(x 1 -x 2 ). (70) 

The two exponents can be rewritten as: — i(p + q) • (xi — x 2 ) — (i/2)(A — A') • (xi — x 2 ) — (i/2)(A + 
A') • (xi + x 2 ). We now make the change of variables: x = xi — x 2 , y = xi + x 2 , as we did after 
Eq. (|35p . The integral over y produces -£/ 3 <5a,a' phis extra terms which vanish in the large L limit. 
Thus, as in the previous section, a, c and g take the values: 
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The integral over x gives the Fourier transform of the correlator F; accordingly, the transition 
probability reduces to: 

EIT^I 2 = A^|/ l (A)| 2 [e p -(q + A)] 2 e-^ +A +P) 2 ^ 

A 

with A defined as in (|4ip . As we see, the structure is minimally modified from that of Eq. (|40p . 
and the answer of [3] for the reduction rate is still obtained. 



at — sin (at) 



(72) 



VIII. COMPUTATION WITH A NOISE CONFINED IN SPACE 

The calculation of the previous Section shows that the reason why g = also for an outgoing 
wave packet, is because a <5a,a' appears, which arises from the integral over space with respect to 
the variable y = xi + X2. This suggests that the problem can be avoided by considering a noise 
which is confined to a finite region of space. We analyze this case here. 

Let us suppose that the correlation function of the noise is: 

E[N(x,t)N(y,s)} = f(t-s)F( X -y)e-^ 2 / e \ (73) 
where I is an appropriate cut off. We start from Eq. (|68p 



w = ^Y.T. n % n %m^]T, (74) 



\2 

k j 



where the temporal part T is given by Eq. (|55p . The two Kronecker deltas coming from Eq. (|6l 
set: k = q + A + p and j = q + A' + p. In this case, Eq. ([70|) is replaced by: 

E[Af*JVy = n 2 7 ( — ) -L f d Xl f ( i X2e ^(p+q)-(xi-x 2 ) e - 4 (A.x 1 -A'.x 2 ) F(xi _ X2 ) e -(x+y) 2 /^ 2 . 
J \m J LP J L s J L 3 

(75) 

As before, we perform the change of variables: x = xi — X2, y = xi + X2- In integrating over the 
new variables, we use the rule: 



2 / 2 



1 r L r+( L - x ) 

dxi I dx 2 f(xi,x 2 ) = - I dx dy[f(x,y) + f(-x,y)]. (76) 

J -(L-x) 



L _L J K ' ' 2 

2 2 



In our case: 



/(xi,x 2 ) = e *(q+A'+p)-x le - l (q+A+p).x 2F(xi _ X2 ) e -(xi+x 2 ) 2 /« 2 (77) 



a'+a 



3 



<-y~ -^F(x)e-y/' =Y[fi(x hyi ), (78) 



= ^ q+ f^+pJ. x ^ (A ,_ A) . y;? ^ v _ y 2 // 2 

1=1 

fi{xi,yi) = e^+ k ^e^ A '- A ^^=i—e- x ^ 4r ce-y^ 12 , (79) 

V47rrc 
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where we have used the Kronecker delta constraints to replace q + p + | (A + A') by | (j + k). 
Thus we arrive at the following expression: 



- 3%- 



n 2 j 



m 



\2 j 3 r L 



m J 8L 6 



x e J(A'-A) tW 



e 



n / dx 



+(L-Xi) 

dt/i 2 cos 



1 



1 -* 2 J^ Ce -yi/l\ 



(80) 



Since e x '?/ 4r c has a cutoff at ~ <C L, x, never approaches L. So we can write (in the large 
L limit): 



cos 



-(j + k)iXi 



dyj ea 



,i(A>-A) jyje -ypl\ 



(81) 



To summarize, substituting Eq. (fHTj) into Eq. ([71]) and noting the Kronecker delta in lif k , the 
effect of having considered a final wave packet in place of a plane wave, and of having confined the 
noise in space, is that the double Kronecker delta <5k,q+p<5j,q+p is replaced by: 



K * = Y,Y, h *^ A ^ A '^+^+p 6 ^+^+pfl(jz I dyj ) 

A A' j'=l ^ L ' 



(A'-A)-y -y^ 2 



(82) 



One can easily check that when £ = oo, the triple integral reduces to #aa'- Kjk then becomes: 
J2a l^(^)| 2 ^k,q+A+p^j,q+A+p ; which implies j = k. The same happens when £ < oo, but h(A) = 
^A,o> ke. when the final state is a plane wave. Thus both a final wave-packet state and a noise 
confined in space are necessary in order to avoid the factor of 2 term. 

Assuming that f < L (thus extending the integration limits to infinity) , and using the Kronecker 
deltas, we can rewrite Eq. (j82|) as follows: 



1 



K jk = /i*(k - q- p)/i(j - q- p)— 7 / dye* 



(83) 



The diagonal elements of Kj k , which multiply the factor of 2 term, give a contribution proportional 
to: 

= El^-q-p)! 2 ^ / d y e ~ y2/l2 « (£f « L ( 84 ) 

Accordingly, the undesired term becomes negligible. Let us consider the remaining terms. The 
overall transition probability can be written as: 



k j 



(85) 
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where 0^ collects all terms which are not crucial for the subsequent evaluation. As in normal 
experimental situations, we assume that the noise is confined in a region which is much bigger 
than the typical width of the wave packet: I 3> A -1 . We can then write: 

E|T fj | 2 ~ G q+p , q+p ^/ l *(k-q-p)^/ l (j-q-p)^3 J dye^-^ e^ 2 ^ . (86) 

k k—j 

In the above expression, we have used the fact that h is peaked where its argument is 0, and we 
have shifted the second sum from j to k — j. The integral, being proportional to e _ ^ k_ ^ 2//16 ^ 2 is 
appreciably different from zero only for |k - j\£ < 1, i.e. |k-j| < t~ x < A. This means that we 
can take: h(j — q — p) = h((k — p — q) — (k — j)) ~ h(k — p — q). So we have: 

EIT^I 2 ~ G q+p , q+p ^|Mk-q- P )| 2 ^3 J dy^ e f(k-J)-y e -y 2 A 2 = 0q+p , q+p , (87) 

k ' k—j 

where we have used the box normalization to continuum normalization replacement 



to write 

; e i(k-j)-y = ( y( y ). (89) 



8L 3 ^ 

k-j 

So this term is not suppressed for £ <C L, provided that we take: L S> rc and L S> £ S> A -1 . All 
these conditions are consistent with typical experimental situations. In this regime, the extra term 
found in [3] is negligible, and the Golden Rule formula used as the basis for the calculations of [1], 
[2] gives the entire answer. 

In terms of continuum state normalization, the point g = that gave the extra contribution from 
the first line of Eq. (|55p is a set of measure zero, and so this point does not give a contribution 
to dT/dt in the large t limit. Since the dependence on I has dropped out, the final answer is 
independent of the assumed spatial cutoff on the noise. So once the box size L — > oo, we can 
then take the limit / — > oo and eliminate the noise spatial cutoff /. We pose the question: If the 
calculation is repeated without a spatial cutoff on the noise, but with the initial and final electron 
wave functions taking the interaction with the noise into account (in analogy with the distorted 
wave Born approximation), will the extra term found in [3] then be suppressed? 



IX. SPONTANEOUS EMISSION FROM THE VACUUM 

The Feynman rules suggest that also the following process is possible: 
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which corresponds to a photon emitted from the vacuum. The analytical expression for such a 
process contains two internal particle's propagators -F12-F21 ^ 0{t2 — — £2) (see Eq. ([13]) ). 

giving a zero contribution. Thus at the non-relativistic level there is no spontaneous photon 
emission from the vacuum. However, such a process is expected not to vanish at the relativistic 
level. 
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Appendix: Noise in the box 

Let us consider a Gaussian noise, with zero mean and correlation function: 

E[N(x,t),N{y, a )] = 6(t-s)F{x-y). (90) 
Writing it in Fourier components: 

N(x,t) = ^3 J dke ikx iV(k,t), (91) 
we easily find the following relation for the correlation function in momentum space: 

E[N(k,t),N(k',s)} = (27r) 3 <5(t-s)<5(k + k')F(k), (92) 

where: 

F(k) = J dxe- ik x F(x) (93) 

is the Fourier transform of the spatial correlator. In placing the noise in a box of size L, we select 
only the Fourier components with the correct boundary conditions, therefore we define: 

+00 



1 T . /?7T \ 



(94) 



j=-oo 
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From Eq. (|92p on can write the correlation function of Nl (j,t): 

E[N L (j,t),N L (j',s)} = L 3 5(t - s )5^yF( T i), (95) 
from which one finds the following correlator for the noise in the box Ni{x,t): 

E[N L (x,t),N L {y,a)] = d(t-s)F L (x-y'), (96) 

with: 

^(x-y) ee ± £ e ^Kx-xO j p ( ^ j) . (97) 

j=-oo 

One can easily prove that in the limit L — > oo, the noise iV^(x, t) as defined in Eq. (|94p . converges 
to N(x,t), and the correlation function Fl(x — y) as defined in Eq. (|97p . converges to F(x — y). 



[1] Q. Fu, Phys. Rev. A56, 1806 (1997). 

[2] S.L. Adler and F.M. Ramazanoglu, J. Phys. A 40, 13395 (2007). 

[3] A. Bassi and D. Diirr, J. Phys. A 42, 485302 (2009). 

[4] S.L. Adler, J. Phys. A 40, 2935 (2007). 

[5] S.L. Adler and A. Bassi, Science 325, 275 (2009). 

[6] G.C. Ghirardi, P. Pearle and A. Rimini, Phys. Rev. A 42, 78 (1990). 

[7] L. Disi, Phys. Rev. A 40, 1165 (1989); 42, 5086 (1990). 

[8] S.L. Adler and A. Bassi, J. Phys. A 40, 15083 (2007). See also references therein. 

[9] G.C. Ghirardi, A. Rimini and T. Weber, Phys. Rev. D 34, 470 (1986). 



